Astronomy & Astrophysics manuscript no. 3706 


February 2, 2008 


(DOI: will be inserted by hand later) 





Number counts in homogeneous and inhomogeneous dark 

energy models 



N. J. Nunes , A. C. da Silva ''^ and N. Aghanim 



o 
o 

(N 

D 
tin 



(N 
> 

o 

\o 

o 
vn 
o 

o 



X 



^ School of Physics and Astronomy, 116 Church Street S.E., University of Minnesota, Minneapolis, Minnesota 
55455, U.S.A. 

^ lAS-CNRS, Batiment 121, Universite Paris Sud, F-91405, Orsay, France 

^ Centro de Astrofiisica da Universidade do Porto, Rua das Estrelas, 4150-764 Porto, Portugal 
February 2, 2008 

Abstract. In the simple case of a constant equation of state, redshift distribution of collapsed structures may 
constrain dark energy models. Different dark energy models having the same energy density today but different 
equations of state give quite different number counts. Moreover, we show that introducing the possibility that 
dark energy collapses with dark matter ( "inhomogeneous" dark energy) significantly complicates the picture. We 
illustrate our results by comparing four dark energy models to the standard A-model. We investigate a model 
with a constant equation of state equal to —0.8, a phantom energy model and two scalar potentials (built out of a 
combination of two exponential terms). Although their equations of state at present are almost indistinguishable 
from a A-model, both scalar potentials undergo quite different evolutions at higher redshifts and give different 
number counts. We show that phantom dark energy induces opposite departures from the A-model as compared 
with the other models considered here. Finally, we find that inhomogeneous dark energy enhances departures 
from the A-model with maximum deviations of about 15% for both number counts and integrated number counts. 
Larger departures from the A-model are obtained for massive structures which are rare objects making it difficult 
to statistically distinguish between models. 
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1. Introduction 

In the present general picture of cosmology, converg- 
' ing evidences suggest that the matter density parame- 
[ ter is low and that the largest fraction of the energy 
density of the universe has an unknown nature leading 
to an accelerating phase. These indicati ons come pr i- 
marily from supernovae la data (e.g . .Riess et al\ l|l998tl : 
iPerlmutter et "dl l|l999|) : lRiess a^J l^2004^ ) and are cor- 
roborated by c osmic microwave background radiation (e.g. 
ISpergel et ai\ tOO^ ) and large scale structure observa- 
tions (e.g lCole etd\ l)2005(l ). alt hough there are differen t 
interpretation s for the data (e.g. iBl anchard et all (|2P0vl); 
I Shanks! l)2004j) '). A cosmological constant can explain the 
acceleration of the universe, however, the disagreement by 
~ 120 orders of magnitudes with predictions from theo- 
retical particle physics has shown the need for resorting 
to an alternative explanation. This is why several theo- 
retical models were recently proposed to explain this dark 
energy in the universe. This new component can be iden- 
tified with a slowly v arying, self - interacting, neutral scalar 
"quint e ssence" field jW etterich' ("1988^ ; 'Ratra and Peebles 
l)l988(l : iFerreira and Jovce (,,1997.) : .Zlatev et aL (,1999) ) 



whic h can b e min i mally coupled, non - minimally coupled 
(e^ UzanI l)l999() : lAmendolal tOO(t: iBaccigalupi et ai\ 
i200.0|)), a phantonTle .g. ICaldwe'iri)l999j) ). a tachvon 
(e.g. Bagla et al! ('20031V) or of purely kinetic nature (e.g. 
lArmendariz-Picon et gj] l)2000(l '). known as K-essence. An 



alternative to dark energy is a fluid with a Chaplygin g as 
typ e of equation of state (e.g. lKamenshchik et a/.l l)200l|) b 
See ISahnil l)2004|) and references therein for a review on 
some of these, and other, models. 



Regardless of its nature, dark energy as a domi- 
nant component, plays a role in the structure forma- 
tion and thus is likely to modify the number of formed 
structures. The evolution of linear perturbations in a 
scalar field like quintessence and the effects on structure 
form ation have already been investigated theoretically 
e.g. IFerreira and Jovcd lll997li : IPerrotta and Baccigalupil 
'2OO2I): lAmendolal 1)2003 )). The effects on the abundance 
of collapsed structures and its evolution with redshift were 
also widely explored and suggested as a to ol to constrain 



the dark energy's nature and evolution (e.g. Haiman et^ 
(2001) : Weller et aL (,2002,) : . Weinberg and Kamionkowski 



2 



Nunes, da Silva & Aghanim: Number counts in homogeneous and inhomogeneous dark energy models 



ll200^. iBattye and Welled tOO^: IWang et al\ l|20a 
iMohJ l)2004| ): lHorellou and Bergd 0bo^ ). 



Recently, numerical simulations including a dark 
energy component were performed by several groups 
to complement the analytical computations and to 
study the effects of dark energy at the structure 
level ( e.g. shape of the d a rk ma tt er halo, ma ss func- 
tionl llLinder and Jenkh^ ll2003l) : iLokas a/.l §003); 
iKlvpin et all ((20n3() : iKuhlen et o.l\ l^2004 ). In such stud- 
ies, the scalar field associated with dark energy is as- 
sumed not to have density fluctuations on scales of 
galaxy clusters or below. If dark energy influences the 
perturb ations on small sc a les as proposed for exain - 
ple bv lArbev et all ll200lh. iBean and Magueilcl 1I2OO2I) 



IPadma ng^hg^^m^Choudhurvf " 1)2002!) or iBaela g/l 

(,2003.) . the collapse of structures as well as their prop- 
ert ies will be affected. 

iMota and van de Bruckl (|20^) have shown that the 
properties of collapsed halos (density contrast, virial ra- 
dius) depend strongly on the shape of the potential, the 
initial conditions, the time evolution of the equation of 
state and on the behaviour of the scalar field in non-linear 
regions. This is what we will refer to as the inhomogeneity 
of the scalar field. More recently. iNunes and Motal l)2004|i 
have investigated how inhomogeneous quintessence mod- 
els have a specific signature even in the linear regime of 
structure formation. They have shown that the time of 
collapse is affected by the inhomogeneity of dark energy 
and they have computed the resulting effect on the linearly 
extrapolated density threshold 5c- Moreover, they exam- 
ined the evolution of matter overdensity as a function of 
time varying equation of s tate in homogeneous an d in- 
homogeneous assumptions. iMaor and Lahavl l)2005f) have 
generalized the formalism to allow for a smooth transition 
between the homogeneous and inhomogeneous scenarios. 
They have concluded that, if only matter virializes, the 
final size of the system is fundamentally distinct from the 
one reached if the full system virializes. 

In the pres ent study, we extend the work of 
iNunes and Motal l|2004) to investigate how the 
quintessence field affects the abundance of collapsed 
halos when the field follows the background evolution 
(homogeneous) and more specifically when it collapses 
with the dark matter (inhomogeneous). We compare the 
two assumptions for models with constant equation of 
state and more general cases of time-varying equation 
of state. To compute the structure abundances and 
their evolution wit h reds hift, we use the canonical 
IPress and SchechteJ l)l974|) formalism. Its theoretical 
expression allows us to account for the effects of inhomo- 
geneous quintessence field through 5c^ the growth factor 
as well as the volume element. Finally, we focus on the 
effects of the different models for structures with masses 
ranging between 10^^ and 10^^ ^-"^Mq. 

This paper is organised as follows. In Section 2 we in- 
troduce the fundamental equations that describe the evo- 
lution of the quintessence field and the collapse of struc- 
ture in the homogeneous and inhomogeneous hypothesis. 



In Section 3 we describe the method used to compute the 
number density of collapsed objects (mass function) in 
both of these scenarios, including the case where the equa- 
tion of state of dark energy is allowed to vary with time. 
We give results and discuss the effects of the normalisa- 
tion of the mass function on the predicted number counts 
in Section 4, and present concluding remarks in Section 5. 

2. Theoretical background 

In a spatially flat Friedmann-Robertson- Walker Universe 
the cosmic dynamics is determined by a background pres- 
sureless fluid (dark and visible matter), radiation and dark 
energy. The governing equations of motion are 



H = 
Pb = 

Pdc = 

with 



Pb + Pdc 

■Pb), 
-Pde) , 



Pdc) , 



Pdc) 



(1) 

(2) 
(3) 

(4) 



Here H ^ a/a is the expansion rate of the Universe, a(t) 
is the scale factor, = SttG and and are the en- 
ergy density and pressure of the background fluid, respec- 
tively. In this work the background is taken to be domi- 
nated by non-relativistic matter, hence, p^ — p^^ oc a~^ . 
If dark energy is a perfect fluid its energy density and pres- 
sure are related by the equation of state pdc — w^Pdc and 
Pdc — f^dcPo/a^^'"'''^''- Alternatively, dark energy can be 
described by a dynamical evolving scalar field rolling down 
its potential V{(t)). In this case, its energy density and pres- 
sure are defined as, p^ — (t? jl -I- V and — (j)^ /2 — V, 
respectively where the relation w — p^/ p^, still holds. The 
equation of motion for the scalar field is. 



dV 



(5) 



We explore models of dark energy that have a signif- 
icant contribution at hi gh redshift unlike in odels such as 
the inv erse power law ilZlatev et all \im^ \ or SUGRA 
models terax and Martini l)l999|) '). We compare the results 
of such models with the cosmological constant model with 
w = — 1 and we assume fide = 0.7 today for all models 
considered in the present study. Moreover, we choose the 
present value of the equation of state, wq, and its running, 
dw I dz(z = 0), to be w ithin the current observational lim- 
its llRiess et all ()2004|) ). 

Figure n shows the evolution of the dark energy den- 
sity and equation of state with redshift for the models we 
consider in this paper. We focus on dark energy models 
for wh ich w = —0.8, w — —1.2 (phantom energv. iCaldwelli 
l|l999|) ') and two cases where the dark energy results from 
a slowly evolving scalar field i n a poten t ial wi th two ex- 
ponential terms (2EXP) ( Barreiro et all \200^ '] 



(6) 
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Fig. 1. Evolution of the dark energy density (top panel) 
and equation of state (bottom panel) with redshift for 
the models considered in this paper: w = — 1 (solid line), 
w = —1.2 (long dashed Hue), w = —0.8 (triple dot-dashed 
line), 2EXPi (a,/3) = (6.2,0.1) (dashed line), and 2EXP2 
(a,/3) = (20.1,0.5) (dot-dashed line). 



where Fde represents the dark energy loss inside the cluster 
and the ratio f/r is, related to the Hubble ratio and the 
evolution of the density contrast through 



1 



r a 
r a 3 1 



If dark energy is homogeneous, then 



Edo ^ I ~ ^ ~ ) iP<^'^c +P'icc) 



(8) 



(9) 



However if dark energy is inhomogeneous and collapses 
with dark matter, Fde = and the equation of motion of 
a scalar field inside the cluster is 



j: dV{(t>c) 



(10) 



Finally, the evolution of the linear density contrast 5l is 
determined by 

^2 

5l = -^HSl + — [p„i5l + (1 + 3w)(5do Pdc + ipdcSw] , 

(11) 

where we have defined the linear density contrast in dark 
energy as 5do = ^Pdo/pdo and 5w is the linear perturbation 
in the equation of state. 



3. Mass function 



We choose as in Nunes & Mota (2004) the pairs (a, /?) = 
(6.2,0.1) (2EXPi) and (a,/3) = (20.1,0.5) (2EXP2). Both 
provide an equation of state at present wq — —0.95. The 
equation of state for 2EXPi approaches zero faster than 
for 2EXP2 (see bottom panel in Fig.^. The two models 
differ also by their contributions at high redshift. As seen 
from top panel of Fig. ^ 2EXPi provides a contribution of 
dark energy that is non negligible at high redshifts (fi^ = 
0.1 at z = 5 for 2EXPi whereas Vt^ = 0.02 for 2EXP2 at 
same redshift). 

In this work we use the spherical collapse model to de- 
scribe the gravitational collapse of an overdense region of 
radius r and density contrast E such that l+(5 = Pm_cl Pn\ = 
(a/r)^, where pnic and p„^ are the energy densities of prcs- 
sureless matter in the cluster and in the background, re- 
spectively. We have considered the possibility that dark 
energy also clusters, i.e. pdcc 7^ Pdc- Following Mota and 
van de Bruck (2004) and Nunes and Mota (2004), we study 
two extreme limits for the evolution of dark energy in the 
overdensity region. First, we assume that dark energy is 
"homogeneous" , i.e. the value of pdc inside the overdensity 
is the same as in the background. Second, dark energy is 
"inhomogeneous" and collapses with dark matter. In gen- 
eral terms, the evolution of dark energy inside a cluster 
can be written as in Mota and van de Bruck (2004) 



Pde, = -3-(/Cdec +Pdec) + Tde , 



(7) 



In hierarchical models, cosmic structures form from the 
gravitational amplification of small initial density pertur- 
bations. The time evolution of the structure abundances is 
determined mainly by the rate at which the perturbations 
grow until they reach the collapse, or virialization. 

An analytical com putation, proposed by 
IPress and Schechteil l)l974j) . gives the comoving number 
density of collapsed dark matter halos of mass M in the 
interval dM at a given redshift of collapse, z by 



dn 
dM 



2 Pmo ^ciz) rflncr(M, z) 



TT M a{M,z) 



dM 



exp 



2a(M,z)2 



(12) 



where pnio is the present matter mean density of the uni- 
verse and (5c (2) is the linearly extrapolated density thresh- 
old above which structures collapse, i.e. 5c{z) = Sl{z = 
Zcoi)- In an Einstein-de Sitter Universe, an overdensity re- 
gion collapses with a linear contrast Sc = 1.686 (see e.g. 
IPadmanabhanI l|l993 )). This canonical value wa s used in 
the seminal paper of lPress and Schechted l|l974|) . The hn- 
early extrapolated density threshold has recently been re- 
computed (Nunes & Mota 2004) in a more general case 
accounting for homogeneous and inhomogeneous dark en- 
ergy. They show significant variations of Sc{z) with red- 
shift and with dark energy models considered. In the 
present study, we compute halo abundances using the re- 
sults on 5c{z) obtained by Nunes & Mota (2004), to which 
we refer the reader for further details. 



4 



Nunes, da Silva & Aghanim: Number counts in homogeneous and inhomogeneous dark energy models 



The quantity a{M,z) = g{z)(7]\j is the hnear theory 
rms density fluctuation in spheres of radius R contain- 
ing the mass M and g{z) = 6l{z)/5]:^{z = 0) is the linear 
growth factor. The smoothing scale R is often specified by 
the mass within the volume defined by the window func- 
tion at the present time (e.g. Peebles 1980). In our analy- 
sis the variance of the smoothed overdensity containing a 
mass M is given by 



fAf = erg 



M 

Ms 



-7/3 



(13) 



where Mg = 6 x lO^'^H.jj^h^^MQ, the mass inside a sphere 
of radius Rg — Sh^^Mpc, and o-g is the variance of the 
overdensity field smoothed on a scale of size i?g . The index 
7 is a function of the mass scale and t he shape parame- 
ter, r , of the matter power spectrum ijViana and Liddl^ 



7 = (0.3r + 0.2) 



(14) 



In our study we use F = 0.167 (Spergel et al. 2003). 

For a fixed erg (power spectrum normalization) the 
predicted number density of dark matter halos given by 
the above formula is uniquely affected by the dark energy 
models through the ratio 5c{z)/ g{z). The underlying as- 
sumption of this approach is that the transfer function 
used in the computation of ctm is that of a cosmological 
constant model. This is a good approximatio n at clus- 
ter sca les for homogeneous dark energy models (|Ma et all 
l|l999(^ ). which remains to be theoretically investigated in 
the inhomogeneous hypothesis. 

For a fixed local halo abundance, the important quan- 
tity to consider is thus = Sc{z)/asg{z). We have ver- 
ified that apart from w = —1.2 (phantom energy) all 
homogeneous models have a ratio Sc/crgg below that of 
the A-model. This means that, at higher redshifts, models 
w = —0.8, 2EXPi and 2EXP2 are expected to give larger 
halo densities (whereas the phantom energy model is ex- 
pected to give lower abundances) when compared to the 
cosmological constant model. For inhomogeneous dark en- 
ergy we find that all models, except the w — —0.8 model, 
have larger Sc/crgg than the cosmological constant model, 
and therefore we expect at high redshift lower hallo den- 
sities for this models (higher abundances in the case of 
w = —0.8) when compared to the A-model. 

To obtain the same halo abundance at z = we scaled 
(Tg according to erg = Sdz = 0)(Ts,A/Sc,Aiz = 0), where the 
index 'a' represents the cosmological constant model, and 
we fix (Tg.A = 0.9. Table^lists the values of cg obtained in 
this way. Note the much larger dispersion of erg between 
models in the inhomogeneous case (third column) as com- 
pared to homogeneous dark energy (second column) where 
models require practically the same erg to reproduce the 
present-day halo abundance of the A-model. 

It is also interesting to see that the 2EXPi model re- 
quires the lowest and the w — —1.2 model the highest 
as normalizations to reproduce the same local abundance 



Table 1. Variance of the overdensity field smoothed on 
8/i~^Mpc used in the normalization of number counts to 
the present-day halo abundance in A-model with erg = 
0.9 for homogeneous (second column) and inhomogeneous 
(third column) dark energy models. 







era 


model 


homogeneous 


inhomogeneous 


2EXPi 


0.897 


0.876 


2EXP2 


0.899 


0.886 


w = -0.8 


0.898 


0.857 


w = -1.0 


0.900 


0.900 


w = -1.2 


0.901 


0.942 



of halos. This means that structures form early, at a slow 
rate, if the universe is dominated by phantom dark energy. 
They form late, at a faster rate, in the 2EXPi model. 

In Fig. 121 we plot the redshift evolution of the mass 
function of objects with mass lO^-* h-^ Mq for both homo- 
geneous (top panel) and inhomogeneous (bottom panel) 
dark energy using Eq. ((T^ . As discussed in the previous 
paragraphs the w = —0.8, 2EXPi and 2EXP2 models give 
larger halo abundances than the w = —I model, whereas 
w = —1.2 gives the lowest densities. We find this same 
qualitative behaviour within the mass range of interest 
for this paper, 10^^ — 10^^ Mq, for both homogeneous 
and inhomogeneous dark energy models. The counts were 
normalized so that at z = all models give the same 
halo abundance as the cosmological constant {w = — 1) 
model gives for erg = 0.9. Note that the overdensity con- 
trast at collapse, Sc, is different for different dark energy 
models (see Fig. 2 in Nunes & Mota 2004) and therefore 
the number density of halos at z = is different if we as- 
sume the same erg for all models. The embedded panels in 
Fig. 121 show this situation (zoomed near z = 0), where CTg 
was set equal to 0.9 for all models. Larger differences are 
found in the case of inhomogeneous dark energy, because 
this is where Sc{z = 0) presents larger deviations from its 
value in a A— model (see Nunes & Mota 2004). 

In what follows we assume the halo number densities 
of models normalized to the present-day halo abundance 
of the A-model. In Section ing we discuss the effects of this 
assumption on our results. 



4. Predicted number counts 

In this paper we investigate the modifications caused by 
a dark energy component on the number of dark matter 
halos. We test for one model with a cosmological con- 
stant (w = —1) and for four quintessence models (de- 
fined in Sectionj^J. The computations arc done in the case 
where dark energy is homogeneous and in the case where 
it may cluster, i.e. inhomogeneous dark energy. Moreover, 
we choose to explore the effects on the integrated num- 
ber of dark matter halos in mass bins [Minf,Msup] illus- 
trating different classes of cosmological structures, namely 
1013 _ 1014, 1014-1015, and lO^^ - lO^^ in units of h~^MQ. 
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Fig. 2. Redshift evolution of the modified Press-Schechter 
mass function at M = 10^"* Mq for homogeneous (top 
panel) and inhomogeneous (bottom panel) dark energy 
models. In the main panels models were normalized to 
reproduce the present-day abundance of dark mater halos 
of the A— model with as = 0.9. In the embedded panels 
models were normalized to the same erg — 0.9. Lines are 
the same as for Fig. ^ 
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Fig. 3. Evolution of the volume element with redshift in 
both homogeneous (top panel) and inhomogeneous (bot- 
tom panel) dark energy scenarios. Lines are the same as 
for Fig. n 

We study the effect of dark energy on the number of 
dark matter halos by computing two quantities. The first 
is the all sky number of halos per unit of redshift, in the 
mass bin 



_dN _ 
dz 



dn / 

47r J Mint 



p dn dV 

dMd^dn'''' 



where the volume element is given by dV/dzdU, — 
r^{z)/H{z), with r{z) — H~^{x)dx. The redshift evo- 
lution of dV/dzdfl for different models of dark energy is 
depicted in Fig. 13 Note that the volume element does not 
depend on the hypothesis of dark energy clustering, it is 
thus the same for both homogeneous and inhomogeneous 
dark energy models. Note also that the phantom model 
is the only one having a volume element larger than the 
volume element in the cosmological constant model. 

The second quantity we compute is the all sky inte- 
grated number counts above a given mass threshold, Mini, 
and up to redshift z: 

(16) 

Our knowledge of both these quantities for galaxy clus- 
ters will improve enormously with upcoming (underway 
or planned) cluster surveys operating at different wave- 
bands. These include the P lanck Surveyor sa tellite and 
South Pole Telescope (SPT) l)Ruhl et al\ ho04)) Sunyaev- 
Zel'dovich surveys, the X MM-Newton s erend ipitous X- 
ray cluster survey (XCS) ijRomer et all (|20^) and the 
recently proposed deep multiband optical Dark Energy 
Survey (DES) designed to probe almost the same sky re- 
gion of SPT. 



4.1. Number counts in mass bins 

Figures 01 [S] and show the number counts, A/"^™ = 
dN/dz, obtained from Eq. IjlSfl (top panels) together with 
the difference of counts of the various dark energy models 
to the A-model, ATV'^''^ = TV"'" - TVj'", (bottom pan- 
els). First let us concentrate on the homogeneous dark 
energy case, i.e. on the left panels of these figures. Below 
redshift unity the volume element has the most impor- 
tant role in the integral of Eq. (|15(l . as it increases by 
orders of magnitude. Above this redshift the volume el- 
ement does not vary much and it is the mass function 
that decreases by orders of magnitude. It decreases faster 
for models of larger Sc/crgg. Therefore, for our particular 
models, we expect to see at low redshifts, in decreasing 
order of number of counts per redshift, the following se- 
quence of models: w = -1.2, w = -1, 2EXP2, 2EXPi 
and w — —0.8 and the reverse order at high redshifts. In 
practice this is only true if the maximum of counts occurs 
at a redshift around or greater than unity, i.e after the 
volume element docs not vary much. We can verify that 
this requirement is only satisfied in the lowest mass bin 
10^^ < Af/(/i~^Mo) < lO^"*, hence the low redshift order 
might differ from the one stated in the two largest mass 
bins. Nonetheless, the high redshift description is accurate 
for all mass bins. More generally, the differences between 
counts can be understood in terms of the relative contribu- 
tion to A/"*^"^ of the integrand in Eq. H15() . The counts thus 
depend not only on the volume element and on the growth 
factor but also on the small differences in as ■ At high red- 
shift, models with smaller 6c/asg ratios imply larger halo 
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Fig. 4. Evolution of number counts (top panels) with 
redshift and differences from the w ~ —1 (cosmological 
constant A) model (bottom panels) for objects with mass 
within the range lO^^ < M/{h-^MQ) < 10^'^. Panels on 
the left show results for homogeneous dark energy whereas 
panels on the right show the same models in the inhomo- 
geneous dark energy scenario. Lines are the same as for 
Fig.ID 



abundances and this effect dominates that of the volume 
element. At low redshifts the differences between counts 
from one model to another depend also on the masses of 
interest. As a matter of fact, the number of low mass struc- 
tures is mostly sensitive to the volume element whereas for 
massive structures the number counts become sensitive to 
the normalisation erg. The left lower panels of the figures 
illustrate quantitatively the relative evolution of the num- 
ber counts. For w — —1.2, the figures depict the expected 
excesses and deficits of counts compared to the cosmolog- 
ical constant model below and above redshift zt (defined 
as A/'^'"(zt) = N'^^'^{zt)). Conversely, we have deficits and 
excesses below and above zt , respectively, for all the other 
models. It is also worth pointing out that larger mass bins 
imply a larger ratio 5c/ cr (see Eq. which makes the 

mass function to dominate at lower redshifts in Eq. p5|l 
and consequently zt to move to smaller values as depicted 
in the figures. 

Now, we turn to the comparison between the inhomo- 
geneous and homogeneous hypotheses (right versus left 
panels in Figs. 0] Eland [HJ. Figure His helpful to under- 
stand how the differences arise. Indeed, when compared to 
the homogeneous case, we verify that the quantity Sc/crsg 
is smaller in the inhomogeneous case for w — —0.8 and 
larger for all the other models. Given this, it is clear that 
for w = —0.8 we see, as expected, smaller deficits and 
larger excesses of counts comparatively to the homoge- 



Fig. 5. Same as Fig. 0] for objects with mass within the 
range lO^-^ < M/Qi-^^Mq) < lO^^. 



10^5 < M/(h-iMo) < 10^6 




Fig. 6. Same as Fig. 0]for objects with mass within the 
range lO^^ < M/[h~^MQ) < lO^^. 



neous case and for all the other models the inverse trend. 
The effects of the inhomogeneous hypothesis are more ev- 
ident for high mass bins where count deficits and excesses 
show larger differences between homogeneous and inhomo- 
geneous dark energy. For example, in Fig. |B| we see that 
the w — —0.8 model has ~ 2 times larger count excess in 
the inhomogeneous case. 
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4.2. Integrated number counts 

The integrated number of coUapsed structures above a 
given mass (Eq. is an important observable quan- 

tity. In this section we present results for the integrated 
number counts of structures with masses above Afinf — 
lQ^^h~^MQ and 10"/i~^Mq. These are displayed in the 
upper panels of Figs. |H1 and respectively (here we omit 
displaying results for M > Mjnf = IO^^H^^Mq because, as 
it will be clear below, integrated counts in this mass range 
can be directly estimated from the curves in Fig. O ■ We 
also calculate the difference in integrated number counts 
with respect to the cosmological constant model shown 
in the lower panels of the same figures. As in the pre- 
vious section, we compare the predicted numbers in the 
homogeneous (left panels of figures) and inhomogeneous 
hypotheses (right panels) for the cosmological models we 
consider in this study. 

For each mass bin, the differences between integrated 
counts result from the combination of effects (that act 
on different mass ranges) discussed in the previous sec- 
tion. In our work we have obviously not performed the 
integration in the mass range in Eq. Iinil aU the way up 
to infinity but only to Mgup = lO^^/i^^M©. In the cases 
M > Mini = 10i3/i-iMo and M > Mnf = IO^/i-^Mq, 
integrated number counts reflect mainly the behaviour of 
curves in the lowest and middle mass bin, respectively. 
This is because N{z, M > Afinf) is dominated by the con- 
tribution of the lower bound of the mass integration range. 
Hence, it is legitimate to concentrate on the lowest mass 
when a qualitative description of the integrated number 
counts is concerned. 

We have seen in the Figs.21[Sland|Hlthat in the homo- 
geneous case, the model with w — —1.2 presents excesses 
with respect to the cosmological constant at low redshift 
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Fig. 8. Integrated number counts up to redshift z (top 
panels) and count differences to the w = —1 (A) model 
(bottom panels) for objects with mass M > IQ^^ H^^Mq. 
Panels on the left show results for homogeneous dark en- 
ergy whereas panels on the right show the same models in 
the inhomogeneous dark energy case. Lines are the same 
as for Fig. 



and deficits at high redshift. All the other models give the 
opposite behaviour as their ^c/csff is always lower than 
that of the w = —1 model. Therefore, for w = —1.2 we 
expect the integrated number counts to be larger than 
for the cosmological constant model until the redshift Zt- 
The remaining models must show the opposite behaviour. 
The difference between the integrated number counts of 
a model with respect to the cosmological constant must 
decrease with redshift for z > zt. Eventually above a red- 
shift Zflat, the integrated number counts should become 
constant with redshift because, as we have noted before, 
the number counts (Eq. decrease exponentially with 

redshift hence contributions above zflat become negligible. 
Note that zflat becomes progressively smaller for struc- 
tures with larger mass limits (Afjnf) because a smaller 
number of these objects form at higher redshifts. This sim- 
ply reflects the hierarchical nature of structure formation 
in models described by Eqs. H12|l - H14|) . Models give quite 
different integrated count differences {N — Nj>^) depend- 
ing on the maximum redshift of integration. For homoge- 
neous dark energy, maximum deviations from the A-model 
are generally obtained near zt of the dominant class of 
structures, whereas for inhomogeneous dark energy maxi- 
mum deviations generally occur at much higher redshifts, 

Z ^ Zfia_f 

In the inhomogeneous scenario, the interpretation of 
the redshift dependence of the integrated number counts is 
similar to that of the homogeneous case. As we have seen, 
the inhomogeneous case yields higher excesses and lower 
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deficits for the w — —0.8 model, hence we expect at high 
redshift, a positive difference in integrated counts com- 
pared to the cosmological constant model. This is quite 
visible in the right lower panels of Figs. |H1 and |51 Through 
a similar argument, we expect at high redshift, a negative 
difference of integrated number counts with respect to the 
cosmological constant for the w — —1.2 model. Moreover, 
because there are hardly any excesses in the 2EXP models, 
we must expect the difference of integrated number counts 
to flatten out at Zt at large negative values for these two 
models. 

4.3. Normalisation 

So far in this work we have normalised the various dark 
energy models such that they reproduce the abundance 
of dark matter halos at redshift zero (i.e. the local abun- 
dance) of a A-model with as ~ 0.9. However, observations 
that are being used to quantify the number of observed 
structures trace baryonic (emitting) gas and not dark mat- 
ter halos directly. For systems where non-gravitational 
physics is important, i.e. the less massive structures, this 
may lead to an important miss-match between dark mat- 
ter halos and emitting structure counts. Therefore, we 
dedicate the rest of this section to discuss the implica- 
tions on our results of dropping the constraint of normal- 
izing models to the same abundance at z = 0. Instead, we 
consider that all models have the same as — 0.9 normal- 
ization, as given by present day observations (e.g. Spergel 
et al. 2003). 

We illustrate the effects of inhomogeneous dark energy 
on one single mass bin. Figure [TUI shows number counts in 
the mass bin 10^'' < M/{h-^MQ) < lO^^ and Fig.lTTlthe 
integrated number counts above mass M > lO^^/i^^M©. 



Both figures consider inhomogeneous dark energy mod- 
els. In the left panels of these figures we are assuming 
that models are normalized to the same halo abundance 
at z = (i.e. same curves as those in the right panels 
of Figs. ElandEl), whereas right panels assume that all 
models have a fixed normalisation, as = 0.9. Models in 
the homogeneous hypothesis have practically the same as 
when they are normalized to the halo abundance at z = 0. 
Therefore we do not expect to observe much difference in 
the halo abundances from one dark energy model to an- 
other. The situation is quite different when dark energy is 
inhomogeneous. The comparison between panels in Fig. 1101 
indicates that fixing as in all models causes much larger 
departures from the A-model than in the case where mod- 
els are normalized to reproduce the same local halo abun- 
dances. At the maximum of A/"*^'", the differences between 
dark energy models come from the fact that for a fixed 
as, the mass function reflects the variations of Sc/g with 
z, which effects dominate those of the volume element in 
Eq. I|15|l for this mass bin. It is interesting to note that the 
structure of the curves can change dramatically depending 
weather we fix the local abundance or as- 

As this work was being completed iManera and Motal 
("2005) made public an analysis similar to the one 
presented here for the particular scenario of coupled 
quintessence. In their work, they have fixed as for the 
two models under study rather than the local halo abun- 
dance. We have verified that indeed, also in the coupled 
quintessence model, the departures from a cosmological 
constant model are larger if one takes the former approach. 



5. Conclusions 

Number counts of collapsed structures are commonly pro- 
posed as a tool to probe dark energy models. In the sim- 
plest case of a constant equation of state, galaxy clus- 
ter number counts may constrain the dark energy, pro- 
vided we have a good knowledge of the cluster physics 
and their redshifts (see for example iMaiumdar and MohJ 
l|2004l) : IWang et all 1^2004) ). In the present study we have 
investigated, using the Press-Shechter mass function, two 
complications to the general picture: First, we have revis- 
ited the modifications introduced by a varying equation 
of state or by including a phantom component; second we 
have explored the effect of an inhomogeneous dark energy 
which collapses with the dark matter. 

More specifically, we have considered a scalar potential 
built out of two exponential terms for which two sets of 
parameters (models 2EXPi and 2EXP2) were explored. 
Although their equations of state at present are almost 
indistinguishable, they undergo quite different evolutions 
at higher redshifts and generally give different results. In 
our analysis we have also included a phantom dark energy 
model (w = —1.2), compatible with present observations. 
An interesting feature of phantom dark energy is that it in- 
duces opposite departures from the A-model as compared 
with the other models considered in this paper. That is, 
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Fig. 10. Redshift dependence of number counts (top pan- 
els) and count differences relative to the A-model (bot- 
tom panels) for objects with mass within the 10^^ < 
M/{h~^MQ) < 10^^, assuming dark energy is inhomoge- 
neous. Panels on the left show results assuming models are 
normalized to the same halo abundance at z = 0, whereas 
panel on the right assume that models have all cts = 0.9 
(which implies different halo abundances at z = 0, see 
text). Lines are the same as in Fig.^ 



we expect an excess of sources in a phantom energy model 
when other models predict a deficit. 

Recently Mota and van de Bruck (2004) proposed that 
dark energy may cluster in forming structures (inhomoge- 
neous dark energy). They have investigated the growth 
and collapse of cosmological structures under the inho- 
mogeneous hypothesis. Going a step forward, we have in- 
vestigated in the present study the implications of this 
hypothesis on the number density of collapsed objects. 
We have found that inhomogeneous dark energy generally 
enhances departures from the A-model. This includes the 
models with a time varying equation of state, which can 
present several times larger departures (from the A-model) 
as compared to the homogeneous case. Yet our results 
indicate that the inhomogeneous dark energy hypothesis 
causes maximum deviations no larger than ^ 15% in mass 
bins with comfortably large numbers of collapsed halos. 
Another interesting feature is that maximum departures 
from the A-model are generally obtained at higher redshift 
for inhomogeneous dark energy than for the homogeneous 
case, which generally show maximum departures near the 
maximum of J\f^™{z). This may be a helpful feature to 
test for the inhomogeneous hypothesis. Larger departures 
from the A-model are also stronger for the more massive 
structures, but these are quite rare objects, which makes 
it difficult to statistically distinguish between models. Our 
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Fig. 11. Integrated number counts up to redshift z (top 
panels) and integrated count differences relative to the 
A-model (bottom panels) for objects with mass above 
10^^/i~^Mo, assuming dark energy is inhomogeneous. 
Panels on the left show results assuming models are nor- 
malized to the same halo abundance at z = 0, whereas 
panels on the right assume that models have all = 0.9. 
Lines are the same as in Fig. ^ 



analysis reveals that the inhomogeneous dark energy hy- 
pothesis has the greatest impact on the 2EXP1 model. 

In this work, we have assumed that the matter trans- 
fer function remains unchanged at cluster scales. We have 
further assumed that models are normalized to reproduce 
the same abundance of dark mater halos at redshift zero. 
In the homogeneous hypothesis all models have practically 
the same erg and there are not much differences in the halo 
abundances from one dark energy model to another. When 
dark energy is inhomogeneous, cts differ by a few percent 
and the departures from the A-model are much larger. It 
is, however, worth noting that the gas physics which rules 
the observed quantities adds a degree of degeneracy. We 
have evaluated the effects of alternatively fixing erg to a 
specific value regardless of the dark energy model. In this 
case we verified that the departures from a A cosmology 
are further enhanced. 

Our results show that constraining dark energy mod- 
els from structure counts is complicated when models 
have time varying equation of state. It becomes an even 
more complicated task when the possibility of inhomo- 
geneous dark energy is taken into account. Therefore in 
order to constrain dark energy models, we need to ex- 
plore as many observable quantities as possible. Our re- 
sults suggest that besides redshift distribution of struc- 
tures, considering structures in mass ranges significantly 
increases the number of observables. Indeed, each theo- 
retical model provides specific predictions for the redshift 
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evolution of number counts and integrated count differ- 
ences in different mass bins. The comparison of such quan- 
tities with observations can be used for testing models 
against observational data. It may further allow to dis- 
tinguish between homogeneous and inhomogeneous dark 
energy models. However, this requires good knowledge of 
the gas physics, redshifts of observed structures and, more 
precisely, a good understanding of the selection function 
of the observations. 
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